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1 Executive Summary 

This report describes the development of a tunable electromechanical Helmholtz 
resonator (EMHR) for engine nacelles using smart materials technology. This effort 
addresses both near-term and long-term goals for tunable electromechanical acoustic liner 
technology for the Quiet Aircraft Technology (QAT) Program. Analytical models, i.e. 
lumped element model (LEM) and transfer matrix (TM) representation of the EMHR, 
have been developed to predict the acoustic behavior of the EMHR. The models have 
been implemented in a MATLAB program and used to compare with measurement 
results. Moreover, the prediction performance of models is further improved with the aid 
of parameter extraction of the piezoelectric backplate. The EMHR has been 
experimentally investigated using standard two-microphone method (TMM). The 
measurement results validated both the LEM and TM models of the EMHR. Good 
agreement between predicted and measured impedance is obtained. Short- and open- 
circuit loads define the limits of the tuning range using resistive and capacitive loads. 
There is approximately a 9% tuning limit under these conditions for the non-optimized 
resonator configuration studied. Inductive shunt loads result in a 3 degree-of-freedom 
(DOF) system and an enhanced tuning range of over 20% that is not restricted by the 
short- and open-circuit limits. Damping coefficient measurements for piezoelectric 
backplates in a vacuum chamber are also perfonned and indicate that the damping is 
dominated by the structural damping losses, such as compliant boundaries, and other 
intrinsic loss mechanisms. Based on models of the EMHR, a Pareto optimization design 
of the EMHR has been perfonned for the EMHR with non-inductive loads. The EMHR 
with non-inductive loads is a 2DOF system with two resonant frequencies. The tuning 
ranges of the two resonant frequencies of the EMHR with non-inductive loads cannot be 
optimized simultaneously; a trade-off (i.e., a Pareto solution) must be reached. The 
Pareto solution provides the information for a designer that shows how design trade-offs 
can be used to satisfy specific design requirements. The optimization design of the 
EMHR with inductive loads aims at optimal tuning of these three resonant frequencies. 
The results indicate that it is possible to keep the acoustic reactance of the resonator close 
to a constant over a given frequency range. An effort to mimic the second layer of the 
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NASA 2D0F liner using a piezoelectric composite diaphragm has been made. The 
optimal acoustic reactance of the second layer of the NASA 2DOF liner is achieved using 
a thin PVDF composite diaphragm, but matching the acoustic resistance requires further 
investigation. Acoustic energy harvesting is achieved by connecting the EMHR to an 
energy reclamation circuit that converts the ac voltage signal across the piezoceramic to a 
conditioned dc signal. Energy harvesting experiment yields 16 mW continuous power for 
an incident SPL of 153 dB. Such a level is sufficient to power a variety of low power 
electronic devices. Finally, technology transfer has been achieved by converting the 
original NASA ZKTL FORTRAN code to a MATLAB code while incorporating the 
models of the EMHR. 

Initial studies indicate that the EMHR is a promising technology that may enable low- 
power, light weight, tunable engine nacelle liners. This technology, however, is very 
immature, and additional developments are required. Recommendations for future work 
include testing of sample EMHR liner designs in NASA Langley’s nonnal incidence 
dual-waveguide and the grazing-incidence flow facility to evaluating both the impedance 
characteristics as well as the energy reclamation abilities. Additional design work is 
required for more complex tuning circuits with greater performance. Poor 
electromechanical coupling limited the electromechanical tuning capabilities of the proof 
of concept EMHR. Different materials than those studies and perhaps novel composite 
material systems may dramatically improve the electromechanical coupling. Such 
improvements are essential to improved mimicking of existing double layer liners. 

2 Objectives 

The objective of this study is to develop electromechanical acoustic liner technology 
for future tunable systems using smart materials technology. Specifically, to achieve this 
objective, the following tasks are required: 

• Develop a physics-based impedance prediction model to convert the selected 
electromechanical inputs into acoustic impedance 

• Experimentally verify the perfonnance of the EMHR and validate the modeling 

• Develop the architecture to enable optimization studies to aid in the design of 
liners to achieve a target impedance range 
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• Design a single cell of an electromechanical liner with the goal achieving a 
NASA-specified impedance range 

• Incorporate the impedance prediction model into the existing NASA-LaRC 
Zwikker-Kosten Transmission Line Code (ZKTL) 

3 Background Information 

Equation Section 1 Originating from the operation of military and commercial 
airplanes, aircraft noise exposure can extend many miles beyond the boundaries of an 
airport. Aircraft noise is potentially dangerous for humans, i.e., aircrew and passengers, 
ground crews and mechanics, and people working or living in the vicinity of the airport. 
With more and more people suffering adverse effects due to aircraft noise, public 
awareness of and concerns for noise emitted from aircraft have increased tremendously. 
As a consequence, worldwide policies and laws have been enacted to restrict the noise 
emission from aircraft. 1 Hence, modern aircraft must reduce their noise emissions to 
satisfy new government regulations. 

2 

Aircraft noise is mainly generated from the airframe and propulsion systems. 
Propulsion noise dominates during take-off (with or without cutback in thrust) and the 
cruise phases of the flight. Airframe noise, on the other hand, is dominant during 
approach. To suppress the propulsion noise, acoustic liners are used to line the wall of 
the turbofan engines. Acoustic liners provide a complex impedance boundary condition 
for noise propagating with the engine duct. The liners can either be passive or active in 
terms of their noise suppression characteristics. A typical passive single degree of 
freedom (SDOF) acoustic liner cell is essentially a Helmholtz resonator that consists of a 
solid backplate, perforated face sheet and honeycomb core. The conventional passive 
SDOF liners are most effective over a narrow frequency range. Multiple DOF liners 
have additional perforated septum sheets and offer a wider suppression bandwidth but 
represent a tradeoff in terms of size, weight, and cost. 4 The greatest limitation of passive 
liners is the constraint of fixed impedance for a given geometry. For a given aircraft 
propulsion system, there may be different optimum liner impedance distributions for the 
differing mean-flow and acoustic source conditions associated with take-off, cut-back, 
and approach conditions. 4 
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Active/adaptive liners have attracted more attention from researchers recently 
because of their promise in suppressing the engine noise under different operating 
conditions. The typical operating conditions of the engine include take-off/cutback, 
cruise, and approach. When operated at different conditions, the characteristics of the 
aircraft noise may be significantly different. Potential ways in which active/adaptive 
liners can modify acoustic liner performance in-situ include changing the geometry 
characteristics of the liner and using bias flow through the liner resistive elements. Little 
et al. tuned their resonator by changing the neck-sectional area. 5 De Bedout et al. 
developed a Helmholtz resonator with a cavity that allows a continuously variable 
volume. 6 The variable volume actuation of the cavity is realized by rotating an internal 
radial wall based on a tuning (control) algorithm. Both methods above mechanically 
modify the geometry of the resonator-like acoustic liner. Generally, there is an obvious 
need for an additional actuator, controller, and power supply. Thus, it is difficult to apply 
these methods to the acoustic liner in an engine. 

Earlier in 1957, Thurston et al. indicated that steady air flow through a orifice will 
change its impedance. 7 In the 1970s, Howe developed a model of the Rayleigh 

o 

conductivity of an aperture, which a high Reynolds number flow passed through. The 
interaction between the incident sound wave and the mean bias flow results in the 
periodic shedding of vorticity. Consequently, a portion of the acoustic energy was 
dissipated into heat. Howe’s work indicates that it is possible to enhance the sound 
absorption using a bias flow through the orifice of the acoustic liner. Hughes and 
Dowling further studied the vortex shedding mechanism for slits and circular 
perforations. 9 ' 10 The perforated liner was experimentally investigated by Sun and Jing. 11 
They found that a bias flow can enhance sound absorption and extend the effective 
bandwidth of the perforated liner. 

In essence, the techniques described above seek to enhance the noise absorption by 
modifying the impedance boundary conditions of the liner system. Furthermore, it is 
worth noting that all adaptive/active systems discussed above, whether by modifying the 
geometry of the acoustic liner or introducing a bias flow to the liner, may require 
actuators, sensors, and controllers. Consequently, such a system is often complex in 
comparison with a passive liner system. Ideally, an acoustic treatment system should be 
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robust and lightweight, have wide noise suppression bandwidth, and have the ability to 
modify the performance of the system in-situ. 

Another type of acoustic treatment using a Herschel-Quincke tube merits some 
discussion here. A Herschel-Quincke (HQ) tube is essentially a hollow side-tube that 
travels along the axis of a main duct (but not necessarily parallel to) and attaches to the 
main-duct at each of the two ends of the tube. " In essence, the HQ tube is a simple 
implementation of the interference principle to attenuate sound. The potential of HQ tube 
to attenuate noise for turbo-fan engine was recently investigated by Smith et al. 13 ’ 14 They 
found that HQ tube can be used to attenuate low-frequency broadband noise, “buzz-saw” 
noise, and blade passage frequency (BPF) tones. The HQ tube can be designed to work 
together with passive liners, which is most effective in attenuating high-frequency noise. 
They also presented an adaptive HQ tube, whose resonant frequency can be tuned by an 
internal throttle-plate flap. Again, additional controller, sensor, and actuator hardware 
are needed. 

A novel method to tune the impedance of the liner system is presented in this project. 
The primary element of this liner is a Helmholtz resonator with the standard rigid backing 
replaced by a compliant piezoelectric composite diaphragm. 1 '’ The electromechanical 
Helmholtz resonator (EMHR) is shown in Figure 1. The acoustic impedance of the 
resonator is adjusted and additional degrees of freedom added via electromechanical 
coupling of the piezoelectric composite diaphragm to a passive electrical shunt network. 

4 Lumped Element Model and Transfer Matrix of the EMHR 
4.1 Description of the LEM of the EMHR 

At low frequencies where, for example, the acoustic wavelength is much larger than 
the largest physical dimension of the device, a distributed system can be lumped into 
idealized discrete circuit elements. 16 Under this approximation, the linearized continuity 
and Euler equations may be replaced by equivalent Kirchoff s laws for volume velocity 
and pressure drop. 

From a lumped element modeling perspective, the EMHR shown in Figure 1 is very 
similar to the compliant isotropic backplate Helmholtz resonator presented in Horowitz et 
al. The primary difference is that the isotropic backplate impedance is replaced by a 
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two-port network representing the electromechanical transduction of the piezoelectric 

jo 

composite backplate. The EMHR is similar in nature to a synthetic jet actuator 
presented in Gallas et al. 19 However, the piezoelectric backplate is actively driven in the 
latter case, while the EMHR is a passive device whose acoustic impedance is adjusted by 
changing the passive loads. 

The equivalent circuit representation for the EMHR is shown in Figure 2, where P 
and P' represent the incident and diaphragm acoustic pressures, respectively. Similarly, 
Q and Q' are incident and diaphragm volumetric flow rates, respectively. In the 
notation below, the first subscript denotes the domain (e.g., “a” for acoustic), while the 
second subscript describes the element (e.g., “D” for diaphragm). Hence, R aN and M aN 
are the acoustic resistance and acoustic mass of the fluid in the neck, respectively. The 
term C aC is the acoustic compliance of the cavity, while C aD and M aD represent the 
short-circuit acoustic compliance and the acoustic mass of the piezoelectric compliant 
diaphragm, respectively. For the diaphragm, M aDrad is the acoustic radiation mass, R aD 
is the acoustic resistance which includes the acoustic radiation resistance and other 
structural damping losses. Finally, C eB is the electrical blocked capacitance of the 
piezoelectric diaphragm (i.e., when the diaphragm does not move so that its volumetric 
flow rate is zero), (j) is the effective acoustic piezoelectric transduction ratio (in Pa/V ), 
Z eL is the electrical load impedance across the piezoelectric backplate, and V is the 


voltage output generated across the piezoelectric backplate under the excitation of an 
incident acoustic wave. From the equivalent circuit representation of the EMHR, shown 
in Figure 2, the acoustic input impedance of the EMHR is given by 

f | j_2rr \ 
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The parameters of the LEM of the EMHR are listed in Table 1. Appendix A presents the 
extraction of the model parameters for the piezoelectric backplate. The analysis of the 
tuning behavior of the EMHR based on the LEM is given in Appendix B. An 
experimental study of the tuning capabilities of the EMHR is presented in Section 8. 


6 


Table 1. LEM parameter estimation for the neck and cavity of the EMHR 17, 18 ' 20-23 



Acoustic 
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c eB 
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C eB = C el , (l - A ) = 5 — All - k 2 ) where £ r is the relative dielectric 

hp 

constant of the piezoelectric material, £ 0 is the permittivity of free space, 
and k 2 = d 2 / C EF C aD is the electroacoustic coupling factor 


4.2 Development of the transfer matrix of the EMHR 

The transfer matrix (also called the transmission matrix or the four-pole parameter 
representation) method has been widely used in acoustic engineering. " The analysis of 
complicated systems is greatly simplified by the use of the transfer matrix method. A 
complicated system can be simply regarded as a “black box” with one input port and one 
output port. Matrix algebra can then be applied to the general treatment of such networks. 
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The interconnection of parallel, series-parallel, and parallel-series network combinations 
can then be handled by simple linear addition or multiplication of the four-pole parameter 
matrix. 26 

The transfer matrix for an EMHR mounted at the end of plane wave tube (PWT) has 
been developed. As shown in Figure 3, the EMHR mounted at the end of the PWT 
composes of the following four components: 

• area contraction 

• area expansion 

• duct element 

• piezoelectric backplate with shunt network 

Hence, it is convenient to develop the transfer matrix for all these elements, and then 
multiply these matrices sequentially to obtain the transfer matrix of the EMHR as 
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Transfer matrix for cavity of the EMHR 

A D <t> 

A D (f) 




Transfer matrix for 
PZT-backplate 


Notations in the equation above, the transfer matrix development details are presented in 
Appendix C. 

5 Optimization of the EMHR 
5.1 Objective of the optimization 

The synthesis of the optimization design problem for the EMHR aims at providing 
the architecture to enable optimization studies to aid in the design of liners to achieve a 
desired impedance range. The optimization design studies of the EMHR follow two 
paths. The first path is to maximize the tuning range of the EMHR with non-inductive 
loads. The EMHR with non-inductive loads is 2DOF and has two resonant frequencies. 
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Both resonant frequencies can be tuned by adjusting the electrical impedance of the non- 
inductive shunts. The difference between the open-circuit resonant frequency and its 
respective short-circuit counterpart is defined as the tuning range. The optimal design of 
the 2DOF EMHR consequently aims to maximize the tuning ranges of both resonant 
frequencies of the EMHR. The second path is to maximize the tuning range of the 
EMHR with inductive loads. The EMHR with inductive loads is 3DOF and has three 
resonant frequencies. 

5.2 Optimizing single tuning range of EMHR with non-inductive loads 

The EMHR with non-inductive loads has two resonant frequency co x and co 2 , and 
thus has two tuning ranges A.a> x and A co 2 . The single objective design optimization 
problem seeks to maximize either tuning range, An^ or A <x > 2 , which are functions of the 

geometric parameters listed in Table 2. To simplify the problem, it is assumed that the 
material properties of the piezoelectric composite diaphragm, i.e. Lead Zirconate Titanate 
(PZT) composite diaphragm, listed in Table 3 are constant. The optimization is 
formulated as follows: 

• Minimize the Objective function: -Au^ or -A <x> 2 

• Design variables: r, t , R , L , R x , h ,R 2 and h s 

• Constraints: 

1) Geometry constraints that impose physical limitations on the design variables. 
The constraints are based on the size of the test apparatus, available commercial 
piezoelectric benders, and size restrictions on the EMHR. 

2) Frequency constraint that confines the first short-circuit resonant frequency f u of 
the EMHR to a particular range (1200 Hz to 1900 Hz) where noise suppression is 
preferred, while also prescribing an upper limit of 3500 Hz for f 2s 

3) Lower bounds (LB) and upper bounds (UB): LB < \r,t,R,L,R { , h p ,R 2 , h s } < UB , 

where LB and UB of each design variable are listed in Table 2. 

Clearly, both the objective function and the frequency constraint are nonlinear; thus the 
optimal design of the EMHR is classified as a constrained nonlinear optimization 
problem. Mathematically, the single objective optimization of the EMHR is 
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Minimize -Aoj 1 (x) or -Aco 2 (x) , where x is the design variable vector 
Subject to 

R-R 2 < 0; R 1 -R 2 < 0; 

-h p - h s + 10^ < 0; h p + h s - 1CT 3 < 0; 

1200-<z> ls /2;r < 0; (y ls /2^--1900<0; o 2j /2^-3500 < 0 
LB i /x i -1 < 0; x i /UB i -1 < 0; 

where {x,.} = i = 1,2...8 . 


Table 2. Design optimization variables of the EMHR (Unit: mm). 


Description 

Symbol 

Lower Bound 
(LB) 

Upper Bound 
(UB) 

Neck radius 

r 

0.5 

3.5 

Neck thickness 

t 

1 

4.5 

Cavity depth 

R 

5 

15 

Cavity radius 

L 

10 

20 

Piezoceramic radius 

R 

1 

25 

Piezoceramic thickness 

K 

0.05 

1 

Shim radius 

r 2 

5 

25 

Shim thickness 

K 

0.05 

1 


The optimization defined above was first implemented using the fmincon command in 
the Matlab Optimization Toolbox and then verified using a Genetic Algorithm provided 
by Matlab Genetic Algorithm Direct Search (GADS) toolbox. One set of the 
optimization results are listed in Table 4. The initial values of the design variables are 
chosen based on the prototype of the EMHR discussed in Liu et al.~ The initial design 
vector satisfies all constraints and thus is a feasible initial condition. Different initial 
design vectors are explored next. The result of the optimization depends strongly on the 
initial values. However, the collective optimal results show that any improvement in one 
tuning range (Aa> x or A a> 2 ) can only occur by compromising the other tuning range. 

Both tuning ranges A<z>, and A a> 2 cannot maximized simultaneously, as shown in Table 4. 
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Table 3. Material properties of piezoelectric composite diaphragm 




Young’s Modules (Gpa) 

63 


Piezoceramic 

Density (kg/m 3 ) 

7700 

Piezoelectric 

(APC 850) 

Relative dielectric constant 

1750 

composite 


Poisson’s ratio 

0.31 

diaphragm 

Shim 

Young’s Modules (GPa) 

110 


(260 half Brass) 

Density (kg/m 3 ) 

8530 


Poisson’s ratio 

0.37 


Table 4. Single objective optimization results of the EMHR (dimension units in mm). 



Initial value 

Optimizing A®, 

Optimizing A®, 

r 

2.39 

3.50 

3.50 

t 

2.18 

1.0 

1.0 

R 

6.35 

5.96 

15.0 

L 

15.2 

12.6 

16.2 


10.1 

21.8 

16.2 

K 

0.12 

0.75 

0.74 

R 2 

12.4 

23.2 

17.2 

K 

0.17 

0.30 

0.23 

(rad) 

211.7 

1085.3 

1.60 

JS[ 

<1 

1.86% 

9.09% 

0.02% 

A®, (rad) 

607.5 

44.4 

2042.6 

A®, / w 2s 

4.14% 

0.2% 

9.29% 


5.3 Pareto optimization of the EMHR with non-inductive loads 

It is impossible to simultaneously maximize both Aco^ and Aa> 2 for the EMHR with 

non-inductive loads. Therefore, a trade-off approach is pursued. Pareto optimization is 
thus explored to optimize both tuning rages at the time to achieve a Pareto solution. The 
Pareto solution (also called a Pareto optimal) is one where any improvement of one 
objective degrades at least one other objective. Three methods have been explored to 
obtain a Pareto solution for multiple-objective optimization of the EMHR: the s - 
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constraints method, 28 the traditional weighted sum method,' 9 and the adaptive weighted 
sum (AWS) method. Mathematically, the a -constraints method is 
Minimize -A(o 2 (x) , where x is the design variable vector 
Subject to 

-Aco x < a; 

R-R 2 < 0; R x -R 2 < 0; 

-h p -h s + 10 4 < 0; h p +h s - 10 3 < 0; , (4) 

1200-<z> ls /2;r < 0; <y ls /2;r-1900 < 0; <z> 2 v / 2zr — 3 5 00 < 0 
LBjx i -1 < 0; xJUB i -1 < 0; 

where {x,} = | r,t,R,L,R v h p ,R 2 ,h s } i = 1,2. ..8 . The tuning range A co 2 is chosen to be the 

primary objective function which is subject to the original constraints and an additional 
constraint limiting the tuning range Aco x . The advantage of the a -constraint method is 

that it is able to achieve the optimal solution even for a non-convex boundary of the 
Pareto front. The problem with this method is that it is generally difficult to choose a 
suitable a . 

The traditional weighted sum method is used to convert the multi-objective 
optimization to a single objective problem by using a weighted sum of the original 
multiple objectives 

Minimize -far, A a?, (x j + (1 - a x j A <x> 2 (V|| , where x is the design variable vector 
Subject to 

R-R 2 < 0; R x -R 2 < 0; 

~h p ~h s + 10 4 <0; h p +h,~ 10 3 <0; 

1200 -co^J 2ft < 0; (o u / 2^-1900 <0; co 2 J 2^-3500 <0 
LBjx,- 1<0; x i /UB i -1 < 0; 

where 0 < a x < 1 is the weighting coefficient. This method is easy to implement, but it 
has two drawbacks. One is that an even distribution of the weighting factors does not 
always achieve an even distribution of the Pareto front. The other is the weighted sum 
method cannot find the Pareto solution on non-convex parts of the Pareto front. 
Furthermore, the traditional weighted sum method can produce non-Pareto points 
occasionally. However, the AWS method is able to produce well-distributed Pareto 
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solutions. For the region unexplored by the traditional weighted sum method, the AWS 
involves changing the weighting factor adaptively rather than by using a priori weight 

30 

selections and by specifying additional inequality constraints. 

The Pareto solution obtained using these three methods are shown in Figure 4. 
Clearly, optimal solutions obtained using the weighted sum method tends to cluster 
together, as mentioned above. A large portion of the Pareto front is obtained using the s- 
constraints method, which is simply ignored by the weighted sum method. This is due to 
the non-convex properties of the Pareto front. Figure 4 also shows the result of the AWS. 
The results indicate that the AWS addresses the drawback of the conventional weighted 
sum method, and produces the Pareto optimum with a good distribution. The result from 
the AWS method matches the result using the s-constraint method very well. 

5.4 Optimization of the EMHR with inductive loads 

The EMHR with inductive loads has 3DOF and thus has three resonant frequencies. 
However, for the EMHR studied, there are two resonant frequencies located within the 
frequency range of interest, as shown in Figure 5. One possible goal of the optimization 
is to find an optimal EMHR in which all resonant frequencies are within the frequency 
range of interest and the frequency shift between the maximum and minimum resonant 
frequency is minimized. Thus, within a wide frequency range, the variation of the 
reactance of the EMHR is relative small. Mathematically, the task is described as 
Minimize A co 
Subject to 

R-R 2 < 0; R l -R 2 < 0; 

~h p ~h s +1CT 4 < 0; h p +h M - 10 3 < 0; 

\200-co L j2ft < 0; (o 3s /2ft -3000 < 0 
LB i /x i -1 < 0; x i /UB i -1 < 0; 

where \x j } = {r,t,R,L,R l ,h p ,R 2 ,h s ,Z eL }i = 1,2. ..8, and the objective function is defined 
as 

A co = max ( co i ) - min ( a > j ) / = 1, 2, 3 (7) 

Note that the design variables include inductive shunt impedance. Moreover, because the 
optimization goal is to limit three resonant frequencies within a given frequency range, 
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the objective function is thus set to be constant, i.e. a given negative real number. Due to 
highly nonlinearity of constraints, Genetic Algorithm provided by Matlab Genetic 
Algorithm Direct Search (GADS) toolbox is thus explored. One set of optimal results is 
listed in Table 5. A comparison of the acoustic impedance between the initial and 
optimal EMHR is shown in Figure 6. It is clear that the reactance of the optimal EMHR 
varies much less than the initial one. 

Table 5. Optimal result for optimization of the EMHR with inductive loads (dimension units in mm) 



initial value 

Optimizing ®,,® 2 ,®, 

r 

2.4 

3.0 

t 

3.2 

1.0 

R 

6.3 

6.5 

L 

16.4 

37.4 


10.1 

13.8 

K 

0.1 

0.1 

r 2 

12.4 

14.8 

K 

0.2 

0.2 

Z ei (mH) 

50 

108.8 

(rad) 

11422.8 

9424.8 

® 2 (rad) 

14903.7 

10386.1 

® 3 (rad) 

24045.8 

12604.1 
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6 Mimicking the NASA 2DOF Liner 


One potential application for an EMHR is to achieve the performance of an ideal 2- 
layer liner while possessing a smaller depth. If achieved, this could potentially result in 
significant cost savings due to weight and size reductions. To explore this possibility, we 
conducted a study to mimic an ideal NASA 2DOF liner design provided by Mr. Jones. 
Specifically, the goal was to replace the second layer of the liner using a piezoelectric 
composite diaphragm. The nominal specific acoustic impedance of NASA 2DOF liner is 
shown in Figure E- 1 . 

The second-layer of a liner can be modeled as a second order system. In an attempt 
to achieve the desired reactance, a curve-fitting method is adopted to obtain the specific 
acoustic mass and compliance that satisfy the NASA DDOF requirements, as shown in 
Figure E-2 and Figure E-3 in Appendix E, which presents how to mimic the NASA 
2DOF liner step by step. It becomes clear that the composite diaphragm should be light 
and compliant to achieve the goal. The PZT backplate, which is used in the EMHR under 
study, thus cannot satisfy such requirement due to its large density and stiffness. 
Moreover, the shim of the PZT backplate is made of half-hard 260 brass which further 
makes it difficult to mimic the NASA 2DOF liner. A possible choice is to use 
Polyvinylidene fluoride (PVDF) instead of the PZT and substitute a light, compliant 
metal shim for brass shim. In comparison with PZT, PVDF is more flexible and has 
lower density but provides relatively weak electromechanical transduction, see Table 6. 
Thus, a trade-off must be reached to design an EMHR, whether to satisfy a given 
impedance requirement or to obtain a satisfactory tuning range. Except for lower density 
and more compliance, some factors must be taken into account to choose the metal shim, 
such as fatigue, operation temperature and cost. One choice is to use a Titanium alloy 
(Ti90/A16/V4); material properties are listed in Table 7. Figure 7 shows a possible 
design using Titanium alloy shim attached with PVDF film, where the damping loss is 
estimated using Eq.(A-9) with f = 0.02 . Note that the actual damping coefficient may 
differ from this and should be determined via measurement, as presented in Section 8.2. 
The results show good agreement between the design and the nominal requirement for 
the acoustic reactance, while there is a discrepancy for the acoustic resistance. The 
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dimensions of this PVDF composite diaphragm are listed in Table 8. With the modeling 
tools in place, improved designs can be obtained in a more systematic study. 


Table 6. Material properties of PVDF and PZT 



Density 

(kg/m 3 ) 

Young’s 

Modulus 

(GPa) 

Poisson’s 

Ratio 

Piezo Strain 
Constant d 3J 
(pC/N) 

Relative 

Dielectric 

Constant 

PVDF 

1780 

8.3 

0.18 

-22 

10-12 

PZT (APC850) 

7700 

63 

0.31 

-175 

1700 


Table 7. Material properties of the Titanium alloy 



Density 

Young’s 

Yield 

Melting Point 



Modulus 

Modulus 



(kg/m 3 ) 

(GPa) 

(MPa) 

(°C) 

Titanium Alloy 
(Ti90/A16/V4) 

4420 

106-114 

1035-1410 

1600-1650 


Table 8. Dimensions of the PVDF diaphragm 



Radius of the 
shim 

Thickness of 
the shim 

Radius of the 
PVDF 

Thickness of the 
PVDF 

PVDF Composite Diaphragm 

10 

0.045 

7 

0.09 
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7 Energy Harvesting Using EMHR 

The acoustic energy harvester is comprised of two components, the EMHR and the 
energy harvester circuit, as shown in Figure 8a. The EMHR transforms incident acoustic 
energy into time-varying electrical energy across the high- impedance piezoceramic. This 
high impedance ac output voltage requires the addition of a power converter to transform 
the power generated into a form that is suitable for an electrical load. Figure 8b shows 
the model of the electromechanical Helmholtz resonator connected to a rectifier bridge, 
large filter capacitor C rect and resistive load i? Load . Maximum power is delivered to the 

load when the magnitude of the EMHR output impedance is matched to the load 
resistance. Shown in Figure 8c is the electromechanical Helmholtz resonator connected 
to a fly-back converter. The converter provides an electronically controllable input 
impedance and transforms the ac voltage output of the Helmholtz resonator into a dc 
voltage delivered to the load. The experimental study of the energy harvesting is 
presented in Section 8. 

8 Experimental Setup and Results 

8.1 Impedance measurement for the EMHR 

The experimental setup implements the standard two-microphone method for 

■j i 

measuring acoustic impedance, shown in Figure 9. The plane wave tube (PWT) 
consists of a 965 mm long, 25.4 mm by 25.4 mm square duct that permits a plane wave 
acoustic field at frequencies up to 6.7 kHz. 

Four Briiel and Kjasr (B&K) type 4138 microphones were used for simultaneous 
acoustic pressure measurements. Two microphones, labeled as Mic. 1 and Mic. 2, were 
flush mounted in a rotating plug to the side of the impedance tube, as shown in Figure 9. 
The rotating plug is used to correct for calibration errors of the microphones by repeating 
the measurements with the microphones in switched positions. The other two 
microphones were used to measure the pressure amplification frequency response 
function in this experiment. One microphone was flush mounted in the sidewall of the 
resonator cavity to measure the cavity pressure. The second was flush mounted to the 
end face of the impedance tube to measure the total acoustic pressure at the entrance of 
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the resonator. This microphone also served as a reference to monitor SPL at the neck of 
the resonator. In the experiment, the SPL is set to be low enough (<100dB) for nonlinear 
behavior of the EMHR to be neglected. 

All microphones are calibrated with a B&K 4228 Pistonphone when connected to a 
B&K PULSE Multi-Analyzer System Type 3560. The PULSE system served as the 
power supply and data acquisition and processing system for the microphones as well as 
the signal source. A pseudo-random waveform generated from the PULSE system was 
fed through a Techron Model 7540 power supply amplifier to drive a BMS H4590P 
compression driver. The driver, which can produce acoustic waves between 200 Hz and 
22 kHz, was connected to one end of the PWT via a transition piece. The pseudo random 
waveform was bandpass filtered from 300 Hz to 6.7 kHz. Meanwhile, a zoom FFT with 
1500 ensemble averages was perfonned on each incoming microphone signal. 

The specimen used in the acoustic impedance measurement using TMM is an EMHR, 
as shown in Figure 1. The piezoelectric backplate is a commercially available 
piezoceramic circular bender disk from APC International, Ltd. The piezoceramic 
material employed in this bender disk is APC850 lead zirconate titanate (PZT). The 
material properties and dimensions of the piezoelectric backplate Helmholtz resonator are 
listed in Table 3 and Table 9, respectively. A variety of loads were used in the 
experiment to investigate the tuning ability of the resonator. Table 10 lists the different 
loads used in the experiment. 


Table 9. Dimensions of the EMHRs (Unit: mm) 




Case I 

Case II 

Neck 

Radius r 

2.42 

2.42 

Length t 

3.16 

3.16 

Cavity 

Radius R 

6.34 

6.34 

Depth L 

16.4 

9.4 

Piezoelectric 

backplate 

piezoceramic 

Radius R, 

10.1 

10.1 

Thickness h 

0.13 

0.13 

shim 

Radius R 2 

12.4 

12.4 

Thickness h s 

0.19 

0.19 
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8.2 Damping coefficient measurement 

It is difficult to accurately model damping of the EMHR. The damping may arise 
from acoustic radiation, thermo-elastic dissipation, compliant boundaries, and other 
intrinsic loss mechanisms. Generally, the present model employed in this study only 
accounts for the acoustic radiation loss. The damping loss measurement thus provides a 
way to check if the acoustic radiation loss represents the damping of the system. The 
schematic of the experimental setup for the damping coefficient measurement is shown in 
Figure 10. The piezoelectric backplate is placed within a vacuum chamber capable of 
producing a low pressure around 10 Torr (1300 Pa). A LK-G32 high accuracy laser 
displacement sensor is used to measure the vibration response of the plate. The 
repeatability of the sensor is 0.05 pm, and the measurement range is 5 mm. A rectangular 
waveform (1 Hz, lOVp-p) generated from the Agilent 33120A Waveform Generator is 
used to excite the piezoelectric plate. This signal is also used to be the external trigger of 
the Tek TDS5104B oscilloscope which acquires data from the laser displacement sensor. 
In addition, 1024 ensemble averages are performed. 


Table 10. Selected loads matrix used in the experiment to tune the EMHR 


Resistive loads 

(O) 

Capacitive loads 
(nF) 

Inductive loads 
(mH) 

nominal 

measured 

nominal 

measured 

nominal 

measured 

Short 

Short 

Short 

Short 

Short 

Short 

200 

199 

10 

10.3 

33 

31.9(70) 

2k 

1.98k 

47 

37.8 

60 

57.2 (130) 

7.5k 

7.44k 

100 

89.7 

100 

101 (65) 

Open 

Open 

Open 

Open 

Open 

Open 


8.3 Parameter extraction of the piezoelectric backplate 

As discussed in Section 4, it is very important for both LEM and TM to accurately 
model the piezoelectric backplate. The parameter extraction provides a way to validate 
the model of the piezoelectric plate. The experimental setup for parameter extraction of 
the piezoelectric backplate is similar to the damping measurement setup described above 
with two exceptions. The piezoelectric backplate is not placed within the vacuum 
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chamber and the transverse displacements along the radius of the backplate are measured 
to reconstruct the deformation of the plate due to the application of the voltage. When 
the defonnation of the plate is measured, the effective acoustic piezoelectric coefficient, 
d a , and the impedance transformer factor, (f> , and the blocked electrical capacitance, 

C eB , of the piezoelectric backplate can then be deduced and compared with the model. 

8.4 Energy harvesting experiment 

The experimental setup for acoustic energy harvesting is similar to the tunable liner 
experimental setup mentioned above, shown in Figure 9, with one exception. The 
electrical tenninals of the EMHR are connected to a power energy harvesting circuit via a 
1 : 1 isolation transfonner which is used to prevent ground loops. 

8.4 Experimental results and discussion 

8.4.1 Acoustic impedance measurement results 

Figure 1 1 shows the acoustic impedance measurement results for the resonator (Case 
I) with a variety of capacitive loads. The EMHR with capacitive loads has two resonant 
frequencies, located where the reactance of the EMHR crosses zero with positive slope. 
There is also an anti-resonant frequency where the reactance of the EMHR crosses the 
zero with negative slope. As the capacitance is increased, the second resonant frequency 
shifts towards the short-circuit case from the open-circuit case, while the first resonant 
frequency barely changes for the EMHR (Case I). This is because of the weak coupling 
between the Helmholtz resonator and the PZT backplate for the tested EMHR. There is 
an approximately 9% capacitive tuning range for Case I resonator under the conditions 
and geometry listed in Table 9. 

Figure 12 shows the measured acoustic impedance for different resistive loads across 
the resonator (Case I). The EMHR with resistive loads has 2DOF. A trend similar to the 
capacitive tuning can be observed for the resistive tuning, in which the tuning range is 
also defined by the short- and open-circuit limits. When the resistance is increased, the 
resonant frequency shifts toward the open-circuit case from the short-circuit case. 
Moreover, resistive loads increase the system damping and thus reduce the amplitude of 
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the impedance peaks. Maximum damping occurs when resistive loads match the 
magnitude of the output resistance of the EMHR. 

Figure 13 shows the measurement results for inductive tuning of the EMHR (Case I). 
Unlike capacitive and resistive tuning, EMHR with inductive loads has 3DOF and thus 
has three resonant frequencies. However, the third resonant frequency of the EMHR 
(Case I) with inductive loads is out of the frequency range of interest and is not shown in 
Figure 13. The third resonant frequency shifts close to the open-circuit case as the 
inductance is increased. Moreover, inductive tuning is not limited by the short- and 
open-circuit cases. The second resonant frequency shifts to lower frequencies away from 
the short-circuit case as the inductance is increased. For EMHR (Case I), the first 
resonant frequency rarely changes due to weak coupling between the piezoelectric 
backplate and the Helmholtz resonator. Furthennore, a real inductor has resistance as 
well. Thus, the peak values of the acoustic resistance and reactance of the EMHR with 
inductive loads has been reduced. 

The experimental results of the EMHR (Case II) with short- and open-circuits are 
shown in Figure 14. It indicates that the first resonant frequency of the EMHR shifts due 
to coupling between the piezoelectric backplate and the Helmholtz resonator is not as 
weakly coupled as in Case I. This is because the depth of the cavity of the EMHR (Case 
II) is reduced. The acoustic compliance of the cavity consequently decreases. Hence, the 
coupling between Helmholtz resonator and piezoelectric composite backplate increases; 
see Appendix B for details. 

The comparison between measurement data and LEM and TM is shown in Figure 15- 
Figure 18 . Clearly, using the measured damping loss of the PZT backplate, both LEM 
and TM predictions match the experimental results pretty well for the short-circuit 
EMHR. The TM provides a better prediction for the acoustic resistance than LEM does. 
This is because the TM includes the viscous loss effect on the wave propagation within 
the cavity. There are some other factors, as presented in Appendix D, which determines 
when the TM coincides with the LEM. 

The predictions for the normalized specific acoustic impedance do not match the 
experimental data very well for the open-circuit case. One possible explanation for the 
observed discrepancy is that, for the short-circuit case, the piezoelectric backplate is 
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electrically shorted and thus (f> and C eB do not affect the acoustic impedance of EMHR, 
as they do for the open-circuit case. Hence, any inaccuracies in either (f> or C eB will 

impact the predicted results. Some factors do impact the piezoelectric backplate model, 

such as the bond layer between the piezoceramic patch and the brass shim (the model 

assumes a negligible bond layer), any asymmetry in the piezoceramic patch geometry, 

18 . 

and imperfect clamped boundary condition. Figure 15 - Figure 18 also indicate that 
modeling damping loss of the backplate as acoustic radiation resistance severely 
underestimate the actual damping loss. 

8.4.2 Damping coefficient measurement results 


Figure 19 shows the damping coefficient measurement results for piezoelectric 
composite backplate in air and in a vacuum chamber. Using logarithm decrement method, 
and if the first two peaks are chosen, the damping coefficient is calculated 

£*= 0-026 ( 8 ) 

for piezoelectric backplate in air, and 

C = 0.024 (9) 

for piezoelectric backplate in a vacuum chamber. However, no measurements are exact, 
and the calculation of the damping factor on first two peaks may lead to errors. Hence, it 
is perhaps more accurate to estimate the damping coefficient by measuring the 
displacement at two different times separated by a given number of periods. For example, 
let x l and x M be the peak displacements corresponding to the times t x and t M , where 

t i+ 1 = t x +iT (/ = 1,2,3...) and T is the period of the system. Moreover, 24 


x i _ e «2</Vi-4- 2 
*< + i 


from which one has the logarithmic decrement 


i x, 




7+1 


from which one has 


in x t = lnxj -8{i -1) . 


( 10 ) 


(ID 


( 12 ) 
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The plot In jc. versus i should has the fonn of a straight line with the slope -8 if the 
measurement is exact. Hence, the accuracy of the estimation for the damping coefficient 
can be improved using Eq. (12). In particular, N peak displacements (x t ,i = 1,2,3...) are 

chosen. In x t and y. = i — 1 are calculated. Then a strait line in the fonn 

z i =ay i +b (13) 

is to be determined to minimize 

N N 

E(ln X, - z t ) 2 = Yj ( ln x t - a y i - b f • ( 14 ) 

i=l i = 1 

Note that a corresponds to —8 in Eq. (12). The damping coefficient is thus calculated 



Twelve peak displacements are chosen from the damping measurement results of the 
piezoelectric plate in air and using method above obtains the damping coefficient of the 
piezoelectric plate in air 

C, =0.015 (16) 

which differs from the estimation in Eq. (8). The illustration of the detennination of 
damping coefficient using the method above and the comparison between numeric 
fittings and measurement results are shown in Figure 20. 

Following similar procedure, one has the damping coefficient of the piezoelectric 
plate in the vacuum chamber 

£*=0.01 (17) 

The illustration of the detennination of damping coefficient in vacuum chamber and the 
comparison between numeric fittings and measurement results are shown in Figure 2 1 . 

8.4.3 Parameter extraction of the piezoelectric backplate 

As discussed above, there is a large discrepancy between the measurement results and 
the predictions of both LEM and TM for open-circuited EMHRs. One possible reason 
for the discrepancy is that any inaccuracies in either (f> or C eB will impact the predicted 

results for open-circuited EMHR. Parameter extraction provides such way to check if 
LEM parameters of the piezoelectric plate are accurate or not. Figure 22 shows the 
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measured displacement of the piezoelectric backplate due to the application of voltages. 
The results indicated that, when applying the voltage on the plate, the measured 
transverse displacement of the plate is larger than one calculated by the model. 
Consequently, the deduced parameters such as d A and (/) and C eB differ from ones from 
the model, as listed in Table 11. Applying the deduced (f> and C eB to the LEM and TM 

results in good predictions of open-circuited EMHRs, as shown in Figure 23. However, 
the reasons for the discrepancy between the predicted and deduced results from the 
measurement data for the LEM parameters of the piezoelectric plate still require further 
investigations. 

Table 11. Comparison between predicted and deduced LEM parameters of the Piezoelectric 
backplate 



d A 

<f> 

c eB 

Predicted 

48.3e-12 

82.4 

45.7e-9 

Deduced 

71.4e-12 

122.0 

41.2e-9 


8.4.4 Energy harvesting experiment results 

To ensure maximum power transfer, it is first necessary to detennine the optimal 
resistive load for the EMHR. A bridge rectifier, filter capacitor and load resistor are 
connected to the Helmholtz resonator output. The sound pressure level at the entrance of 
the EMHR is held constant at 140 dB, while the load resistor is swept from 1 kQ to 100 
kQ . The incident acoustic power is measured using the two microphone method. The 
electrical power is measured by finding the voltage across the load resistor and the 
current flowing through it. Figure 24a shows the output power of the EMHR versus load 
resistance. The experimental results show that the optimal load for the EMHR is — 18 
kQ . Based on this result, the input impedance of the flyback circuit was set to 1 8 kQ . 
The flyback converter was also compared to the direct charging method, where the 
EMHR is connected to a rectifier bridge and large capacitor with the rechargeable battery 
placed in parallel with the capacitor. The power flowing into the battery is measured and 
this serves as a baseline for comparison. Figure 24b presents comparisons between the 
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flyback converter, the direct charging method. The data indicate that for SPL above 140 
dB the flyback converter circuit is superior to the counterpart method. Specifically, at an 
incident SPL of 153 dB the flyback converter harvests 139% more power than the direct 
charging method. 

9 Conclusions and Recommendations 

The EMHR is a promising technology that may enable low-power, light weight, 
tunable engine nacelle liners. Lumped element model (LEM) and transfer matrix (TM) 
of the EMHR have been developed to predict the acoustic behavior of the EMHR. The 
models have been implemented in a MATLAB code and incorporated into the NASA 
ZKTL code, which is converted to a MATLAB code from the original FORTRAN code. 
In this study, LEM and TM match each other pretty well within the frequency range of 
interest. The difference in the results is due to whether the viscous loss is taken into 
account when modeling the cavity of the EMHR. Moreover, the LEM coincides with the 
TM when some conditions are satisfied, which are listed in Appendix D. The analysis of 
the tuning behavior of the EMHR based on the LEM validates the experiment 
observation. 

An experimental investigation of the EMHR has been implemented using the 
standard two-microphone method (TMM). The measurement results verify both LEM 
and TM of the EMHR. Good agreement between predicted and measured impedance was 
obtained. Short- and open-circuit loads define the limits of the tuning range using 
resistive and capacitive loads. There is approximately a 9% non-inductive tuning limit 
for the second resonant frequency under these conditions for the non-optimized EMHR 
configuration studied. Inductive shunt loads result in a 3DOF system and an enhanced 
tuning range of over 20% that is not restricted by the short- and open-circuit limits. The 
damping coefficient measurement for the piezoelectric backplate in the vacuum chamber 
is also performed. The results show that the acoustic radiation damping loss is relatively 
small. Moreover, parameter extraction of the piezoelectric backplate is 
implemented to improve LEM and TM prediction performance, i.e. when EMHR is 
not short-circuited. 
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Based on models of the EMHR, the Pareto optimization design of the EMHR has 
been performed for the EMHR with non-inductive loads. The EMHR with non-inductive 
loads is 2DOF and have two resonant frequencies. Either of them can be optimally tuned 
by approximately 9%. However, the results show that it is impossible to maximize both 
tuning ranges simultaneously. The improvement of one tuning range degrades the other. 
Consequently, a trade-off must be reached. In other words, a generally accepted Pareto 
solution should be achieved. Three methods are explored to obtain the Pareto optimal set 
of the bi-objective optimization design of the EMHR: the s -constraint method, the 
traditional weighted sum method, and the adaptive weighted sum method. Both the s - 
constraint method and the adaptive weighted sum method obtain the same Pareto front 
for the optimization problem, while the weighted sum method clusters the Pareto 
solutions around the two single-objective design points of the Pareto front. The Pareto 
solution provides the information for a designer that shows how design trade-offs can be 
used to satisfy specific design requirements. The EMHR with inductive-loads is 3DOF 
and has three resonant frequencies. The optimization design of the EMHR with inductive 
loads aims at optimal tuning of these three resonant frequencies, i.e. to constrain three 
resonant frequencies within a given range. 

An effort to mimic the second layer of the NASA 2DOF liner using piezoelectric 
composite diaphragm has been made. The optimal acoustic reactance of the second layer 
of the NASA 2DOF liner is repeated using a thin PVDF composite diaphragm while the 
acoustic resistance of the mimic requires further investigation. The PVDF composite 
diaphragm is more flexible than the PZT composite diaphragm but possesses relatively 
weaker electromechanical transduction. Hence, the EMHR with a PVDF composite 
diaphragm has a small tuning range. Clearly, this is an area in which novel composite 
structures may be able to achieve improved perfonnance whether it be in terms of 
increasing damping or tuning range. 

Acoustic energy harvesting is achieved by connecting the EMHR to an energy 
reclamation circuit that converts the ac voltage signal across the piezoceramic to a 
conditioned dc signal. Energy harvesting experiments yield 16 mW continuous power 
for an incident SPL of 153 dB. Such levels are sufficient to power a variety of low power 
electronic devices. 
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Future work should include the testing of sample EMHR liner designs in NASA 
Langley’s dual-incidence and grazing-incidence flow facilities to evaluating both the 
impedance characteristics as well as the energy reclamation abilities. Inductive tuning 
offers the capability of an increased tuning range, but further work is needed regarding 
the design of such liners. Specifically, system goals are required to drive the formulation 
of the optimal design problem. Poor electromechanical coupling limits the 
electromechanical tuning capabilities of the EMHR. Novel composite material systems, 
like an interdigitated piezoelectric composite may dramatically improve the 
electromechanical coupling. Such improvements are essential to mimic both resistance 
and reactance of existing double layer liners. Furthermore, enhanced electromechanical 
tuning will also enable the use of more complex passive ladder circuits for impedance 
spectra shaping. 
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Figure 1. Schematic of the EMHR 



Figure 2. Equivalent circuit representation of the EMFiR 
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Figure 3. Illustration of the modeling EMHR using the transfer matrix method, where each 
subsystem is denoted as: 1-area contraction; 2-duct element; 3-area expansion; 4-duct element; 5- 
piezoelectric backplate. 
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Figure 4. The Pareto front obtained via the £-constraint, traditional weighted sum, and adaptive 
weighted sum methods. 
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Figure 5. Three resonant frequencies of the EMHR with inductive as function of the inductive loads 
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Figure 6. Comparison between initial and optimal acoustic impedance of the EMHR with inductive 
loads 
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Figure 7. Mimic NASA 2DOF liner 
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(b) 


1:1 



(c) 

Figure 8. (a) Schematic of energy harvesting using EMHR. (b) Converter type I . (c) Fly-back 
converter 
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Figure 9. Schematic of the acoustic impedance measurement for the EMHR. 
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Figure 11. Experimental results for the normalized specific acoustic impedance of the EMHR (Case 
I) as function of the capacitive loads. 
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Figure 12. Experimental results for the normalized specific acoustic impedance of the EMHR (Case 
I) as function of the resistive loads. 
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Figure 13. Experimental results for the normalized specific acoustic impedance of the EMHR (Case 
I) as function of the inductive loads. 
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Figure 14. Experimental results of the normalized acoustic impedance of the EMHR (Case II) for the 
short- and open-circuit. 


42 




1500 2000 2500 3000 

Freq.[Hz] 


Figure 15. Comparison LEM, TR and measurement results for a short- and open-circuited EMHR 
(CASE I), the damping loss of the backpiate is determined using logarithm decrement method 

( 4 " = 0 . 015 ) 
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Figure 16. Comparison LEM, TR and measurement results for a short- and open-circuited EMHR 
(CASE I), the damping loss of the backplate is assumed to be acoustic radiation resistance 
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Figure 17. Comparison LEM, TR and measurement results for a short- and open-circuited EMHR 
(CASE II), the damping loss of the backpiate is determined using logarithm decrement method 

(4" = 0.015). 
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Figure 18. Comparison LEM, TR and measurement results for a short- and open-circuited EMHR 
(CASE II), the damping loss of the backplate is assumed to be acoustic radiation resistance. 
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Figure 19. Damping coefficient measurement for piezoelectric composite backplate (Case I) (a) in 
air (b) in vacuum chamber 


47 


0.4 


0.2 


0 


- 0.2 

N 

x- 

“ -0.4 


- 0.6 


- 0.8 


-1 L 
0 


Data 

z.=-0.09439y.+0.1 061 




-L 


L 


2 4 6 8 10 12 


i 


(a) 



Time [s] x 10 ‘ 3 


(b) 

Figure 20. (a) Curve fitting the measurement data (in air) using a 2 nd -order system, fitting 1- 
£ = 0.026 , fitting 2 -£ = 0.015 . (b) Determination of damping coefficient of the PZT plate in air 
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Figure 21. (a) Curve fitting the measurement data (in air) using a 2“ d -order system, fitting 1- 
£ = 0.024 , fitting 2 -£ = 0.01 . (b) Determination of damping coefficient of the PZT pfate in air 
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Figure 22. Measured transverse displacement of the piezoelectric backplate due to the application of 
various voltages. 
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Figure 23. Predictions of the LEM and TM for short- and open-circuited EMHRs. (a) EMHR (Case 
I). (b) EMHR (Case II). 
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Figure 24. (a) Power delivered by the Helmholtz resonator to a resistive load, (b) Output power 
versus incident pressure for the flyback circuit and direct charging method. 
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Appendix A. Parameters Estimation for LEM of the Piezoelectric Composite 
Diaphragm 

Equation Chapter (Next) Section 1 

The extraction of the model parameters for the piezoelectric backplate is more 
complex due to composite plate mechanics. The piezoelectric backplate consists of an 
axisymmetric piezoceramic of radius R t and thickness h p bonded in the center of a metal 


shim of radius R 2 and thickness h s . Up to and just beyond the first resonant mode, the 
one-dimensional piezoelectric electroacoustic coupling is given by 18 


AV 


q 




(A-l) 


where AV is the volume displacement of the piezoelectric backplate due to the 
application of the pressure P and voltage V . Additionally, q is the charge stored on the 

piezoelectric electrodes, d a is the effective acoustic piezoelectric coefficient, and C eF is 
the electrical free capacitance of the piezoelectric material. The volume displacement is 
calculated by integrating the transverse displacement, w(r ) , over the whole plate 


r 2 

AV = | 2?rrw(rjdr . 
o 


(A-2) 


Thus, the short-circuit acoustic compliance of the backplate, C aD , is determined by 


C aD = 


AV 


“2 

J 2k rw ( r yir 


v=o 


v=o 


(A-3) 


Similarly, the effective acoustic piezoelectric coefficient, d A , is detennined by 
application of voltage to the free piezoelectric plate 


d 


A 


AV 

V 


lp=o 



(A-4) 


The impedance transformer factor, (j ) , is 


|p=o 
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(j) = -^A. 

c 

^ nPt 


(A- 5) 


The effective acoustic mass of the piezoelectric backplate, M aD , is found by equating the 

1 8 

total distributed kinetic energy stored in the velocity of the plate to a lumped mass as 

r 2 


M aD = 


2 n 
AV 


J Pa ( h A')) 2 rdr 


(A-6) 


v=o 


where p 4 is the area density of the piezoelectric backplate. The blocked electrical 
capacitance, C eB , is related to the free electrical capacitance of the piezoelectric 
backplate, C EF , as 


C„ = C BF ( i-S) = SM*L(i-S), 


(A-7) 


where s r is the relative dielectric constant of the piezoelectric material, is the 
pennittivity of free space, is the radius of the piezoceramic, and h is the thickness of 
the piezoceramic, and k 1 is the electroacoustic coupling factor and given by 


k 2 = 


d] 


C C 

'-'EF'^aD 


(A-8) 


The acoustic resistance, R aD , of the PZT backplate models acoustic resistance and 
structural damping in the backplate. The damping may arise from thermo-elastic 
dissipation, compliant boundaries, and other intrinsic loss mechanisms. The acoustic 
resistance is given by 

[m +m 

R aD =2C °° (A-9) 

V ^ aD 

where C, is an experimentally determined damping factor determined using, for example, 

24 

the logarithmic decrement method. 
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Appendix B. Capacitive and Inductive-Tuning Strategy of the EMHR 

Equation Chapter 2 Section 1 Equation Section 2 

Capacitive tuning of the EMHR 

The EMHR with capacitive shunts is analogous to a 2DOF system with two resonant 
frequencies, as shown in Figure B-l. To simplify the problem without loss of generality, 
it is assumed that the EMHR with capacitive shunts has negligible damping. 

Referring to Figure B-l a, by defining 

(B-l) 


and 




(c^+cjc^c, 


aC 


(c, B +c aL )c aD +(c„+ c, L + c aD ) c aC 


(B-2) 


the coupling coefficients, which defines the ratio of the oscillating energy stored in the 
coupling elements to that stored in the total capacitance for each loop, are given by 34 


M/c* c, 
(\Q,dtf/c, C .c 

(la^)Vfic c, 

(Ja*) 7 c > C * c ’ 


(B-3) 


(B-4) 


where Q t is volume velocity, as shown in Figure B-l . Hence, the impedance of loop 1 is 


Z L\ ~ R aN + J 
= R aN + J 


6)M aN - 
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coC 


aC 
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^ aN 

V C aC 


CO 






CO 


aC y 


(B-5) 


= z 
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aN 


wLl 


J L\ 


+ j 


CO co n 
V ® J 


~ Z wL\ { A Ll + J^Ll} 


where the quantity 

<% = 1 (B-6) 

is the resonant frequency of the loop 1 . Next, 
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(B-7) 

is the dissipation factor, which is the ratio of the power dissipated to the power stored. 
The quantity 

=®/au-(°Li/ (0 (B-8) 

is the tuning factor which measures the deviation of the operating frequency of loop 1 
from its resonant frequency. Finally, The quantity 

z*. = jKJeTc (B-9) 


is a weighting factor by which the magnitude of the impedance of a different system 
which has the same dissipation factor and the resonant frequency differs from each other. 
Similarly, the impedance of the second loop is 


D 

y y I v aDrad , • 

^12 ~ Zj wL2 | y ^ J 

^ wL2 


CO 


CD, 


L 2 


CO , 


L2 


CO 


Z W L2 (^L2 + J^L2 ) ’ 


(B-10) 


where the definitions of the quantities co L2 , Z wL2 , A L2 and Q i2 are similar as their 
counterparts in loop 1 

/V( M .» +M .iw)C 2 , (B-ll) 

A i2 =^/V(M„ d +M^)/C 2 , (B-12) 

Q. l2 = co/ (o 2 - o) 2 l co , (B-13) 


and 

Z, I2 = d M ^ +M ^w)/C. (B-14) 


Next, for the coupled system 


and 


P = Z l Q 1 + 


jcoC aC 


q 2 . 


(B-15) 


which leads to 


0 - Z 2 Q 2 +—~ — Q x , 
J«oC aC 


Z =^ = Z - — 

^ alN ^ ^ L\ rj 

^L2 


( , V 


jcoC i 
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(B- 16) 


(B-17) 
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From Eqs. (B-5), (B-10) and (B-17), one has 


7 1(1)1 

-^ = A II +jQ il : : 

Z w l\ Z wL] Z wL2 \jooC aC j A L2 +jCl L2 


If the system is undamped or weak damped, 


^£1 ^£2 ^ • 


Thus, Eq. (B-18) is rewritten as 


(B-18) 


(B-19) 


= 7^i — 


Z W L\Z„L2 U<°C aC ) J Q L2 


i co co IA i co Ll co L2 ^2 

— / i ; 1C 


(0 I A CO J CO 


CO CO, 


co l2 CO 


where S=C,CjCl c = k x k 2 is the coupling factor of the system. At the resonant 


frequency, Eq. (B-20) equals zero, and the solution for the resonant frequencies is 


2 

®1,2 — 


(®li +co L1 ) + J((o u +®z, 2 ) 4^1 k ^((o Ll (o L2 ) 


(B-21) 


Furthermore, Eqs. (B-6) and (B-l 1) yields 


<X) L \ + CO L2 


M M C aC {M, d +M m )C 2 


c + c + c 

^aB ^ ^ aL ^ ^ aD 


M aN M aD + M aDrad _ C aC ( M aD + M aDmd ) ( C aB + C aL ) C a 


(l + oc^ + - 


C + C + C 

^aB ~ ^ aL ^ ^ aD 


(B-22) 


M, n C x ' (M, iD +M m )(C„ i+ C, iL )C^ 

= ( 1 + CC ) (Off B + C Off 

where a = M aN /[M aD +M aDmd ) is the mass ratio between the neck and piezoelectric 
backplate. In addition, 


0J HR 


M aN C aC 


(B-23) 
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is the resonant frequency of the Helmholtz resonator with solid wall instead of a 
piezoelectric backplate, and 


(o,, 


C f n + C j + C n 

aEB aL aD 


(B-24) 


^{M aD +M aDrad )(C aB + C aL )C aD 

is the resonant frequency of the piezoelectric backplate. Similarly, 

(B-25) 

• If the system is weakly coupled, k — > 0 , substitution of Eqs. (B-22) and (B-25) into 
(B-21) results in 


2 2 / \ 2 2 

o Ll — co L1 =yl — a) co UD — co r 


HR 


CO\2 


( 1 + cc ) co HR + co D i a/ ( 1 cc ) CO HR CO ( 
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• If the mass ratio between the neck and the piezoelectric backplate is very small, 
cc — > 0 , 


® 1,2 


CO hr + (0~ D ± 


®HR 


n2 


(B-27) 


co. 


HR 


(Or 


In other words, the resonant frequencies of a lightly damped EMHR with capacitive 
shunts possessing weak coupling are approximately the resonant frequency of the solid- 
walled Helmholtz resonator and the piezoelectric backplate. As indicated by Eq.(B-24) , 
the resonant frequency of the piezoelectric backplate is adjusted with the change of the 
capacitive loads 

’ ' >oo 


o Open-circuit , C aL = 0 or Z aL = l/.vC 

I 

V 0 !) Ln = 


c + c 

^ aEB ^ ^ aD 


(M aD + M aDr ad ) ^ aEB ^ a, 

o Short-circuit, C aL =oo or Z aL = 1 /sC aL — > 0 


(B-28) 


\co n , = 

L Jshort 


1 


[M aD + H aDmd ) C aD 


(B-29) 
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o Capacitive load, C aL = <f~ jC eL 



(B-30) 


Clearly, 



(B-31) 


or by the use of Eq. (B-27) 


\co 7 ] u < [ 1 . . <\co 7 1 

L ^ J snort L ^ J capacitive L ^ Jopen 


I capacitive 


(B-32) 


Eq. (B-32) demonstrates how the tunable electromechanical Helmholtz resonator works 
with different capacitive shunts. The short and open loads define the limits of tuning 
using capacitive loads. The second resonant frequency shifts toward the short case when 
the capacitance increases. It is worthy to note that the aforementioned analysis focuses on 
the capacitive tuning behavior of the EMHR with weak coupling. In fact, when such 
coupling is not weak, the first resonant frequency of the system can be tuned in-situ as 
well, as indicated by Eq.(B-21). 

Inductive tuning of the EMHR 

As shown in Figure B-lc, The EMHR with inductive loads is a 3DOF system. 
Similarly, to simplify the problem without loss of generality, it is assumed that the 
EMHR with inductive shunts has negligible damping. The effective impedance of loop 3 
is thus given by 



(B-33) 


which can be further approximated as 


Z L3 □ j(0(f~M eL 


(B-34) 


with the assumption 


® > 2 C eB M eL « 1 , 


(B-35) 


Or 


co Q 


1 


(B-36) 
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where co E is the resonant frequency of the loop 3. Eq. (B-36) is satisfied for the EMHR 
and the frequency range investigated in this project. Physically, loop 3 effectively adds a 
mass (f> 2 M eL to loop 2. Similarly, the resonant frequencies of the weakly coupled EMHR 
are given by 


co, 


1,2 


CO hr + CO D + 


®HR D 


CO 


HR 


(Or 


(B-37) 


where co HR is given by Eq. (B-23), while co D is 


co r 



1 


(B-38) 


Clearly, co D < [ru^,] in Eq. (B-29). With increasing inductive loads, co D moves further 
away from the short-circuit resonant frequency . 
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Z aL = <t > 2 / 'j®C eL 
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Z aL = j(0(jrM eL 


(a) 


(b) 


(c) 


Figure B-l. EMHR with passive electrical loads is analogous to a 2DOF system for (a) 
capacitive and (b) resistive loads and a 3DOF system with an (c) inductive 

load. 
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Appendix C. Development of Transfer Matrix of the EMHR 

Equation Section (Next) 

Area contraction 

The area contraction is shown in Figure C-la. Due to abrupt change of the area, even 
with the planar incident wave, evanescent higher order acoustic modes will be excited in 
the vicinity of the discontinuity, thus the acoustic field near the discontinuity will have a 
transverse component. Such effect can be taken into account in terms of plane wave 
variables by means of an additional mass and resistance. ’ Following Karal and Ingard, 
the relation between the acoustic variables for the area contraction is given by 

P s =P t +(p a c„(JA t )U,A,, (C-l) 

Po C 0^5 A 5 — Po C 0 U 4 A 4 (C"2) 


where P 4 and P 5 are the plane wave pressure component in before and after the 


discontinuity. Eqs. (C-l) and (C-2) represent the continuity of the pressure and the 
volume velocity across the discontinuity, respectively. The matrix representation for the 
area contraction is thus given by, 


Ps 


"1 

c " 

Pa 

_ p 0 cjj 5 _ 


0 


_ P 0 c 0 U 4 _ 


(C-3) 


where A, = nr 1 is the cross-sectional area of the neck of the Flelmholtz resonator, 


A 5 = Dl ihe is the cross-sectional area of the tube, and c is the normalized specific 


acoustic impedance, which is given by 
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C=- 


1 


Po c o 


V 2 m 


CO + 


A 4 

2 nc , 


+ 


o y 


j%kr 
3 n 


(C-4) 


Thus, the transfer matrix for the area contraction is given by 

"1 C 


T = 

1 AC 


0 a 4 /a 5 


(C-5) 


Area expansion 

Similarly, the relationship between the acoustic variables for the area expansion, 

•3 /T 

shown in Figure C-lb, is given by 

P 3 =P 1 + (p,c a (JA 1 )U 2 A ! , (C-6) 
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Po C 0 ^ 3^3 Po C J^2 J ^2 ’ 


(C-7) 


or in transfer matrix form, 


p 3 


"1 

P 

P q cJJt, _ 


0 4/4 

_Po C 0 U 2_ 


(C-8) 


where A 3 = nr 2 is the cross-sectional area of the neck of Helmholtz resonator, A 2 = nR 2 


is the cross-sectional area of the cavity, and t^ 3 is the normalized specific acoustic 
impedance, which is given by 


£>=- 


Po c 0 


V2 m 


CO + 


P 0 cd 2 A 3 


2nc, 


+ 


0 J 


Mr 

3 n 


Thus, the transfer matrix for the area contraction is given by 


(C-9) 


1 (AMK, 

0 A, /A, 


(C-10) 


Duct element 

Finding the acoustic velocity and pressure fields due to sound wave propagation in a 
long, round duct is a classic acoustics problem. Tij deman reviewed this topic and 
showed that the solution obtained by Zwikker and Kosten has the widest range of 
validity. ’ For plane-wave propagation in a duct with radius R , as shown in Figure C- 


lc, the pressure and axial velocity component are given by 
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P = 


U-cA 


Wo 

V Y . 


I (Ae- rkx + Be rkx ) 


Y 


J o(f 2s ) 


( a -Ykx d TAx \ 

lAe -Be j. 


(C-ll) 


(C-12) 


where T is the complex propagation coefficient, k = c o/c 0 is the wave number, // = r/ R 

is dimensionless radius, 5 =R 3 Joop/ju is the shear wave number or Stokes number, y is 
the ratio of specific heats, and p is dynamic viscosity coefficient of the air. The 
complex propagation coefficient is given by 


r = 



(C-13) 
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where 


n = 


1 + 


r- 


-1 J 2 (j m <rs) 




(C-14) 


r Jo{f 2(j s) 

where a = ^J/iC p /A = VPr , Pr is Prandtl number and ./ ( ) is the i th order Bessel 
function of first kind. 

From Eq. (C-12), the mean velocity over cross-section of the duct is thus obtained by 

1 R 

— t- f U 2 7i r dr 

rR 2 J 


U = 


nR- 


t (Ae- rla -Be rkx ) 


I 


JT 


7lR 2 

jTc„ J_1 (i 1,2 d 

r A (/'"■') 

= ^lLi A e- rkx -Be Tkx ) 
y Tn x ' 


y 


J 0 {j m ns) 


Inrdr 


(C-15) 


/ a —Tkx r) 

lAe -Be I 


where two identities of the Bessel function are used 

In 
z 


J n-x{ Z ) +J nA Z )=— J n{ Z )' 

Z 

a 

J zJ 0 (z'jdz = aJ l (a) . 


(C-16) 

(C-17) 


Thus, from Eqs. (C- 11) and (C-15), one can deduce the relationship between ( P 2 ,U 2 ) 
and , illustrated in Figure C-lc, as 


p 2 


_Po c J~y 2 _ 



cosh(rkE) sinh(rkE)/G 
Gsinh(rkL) cosh(rkL) 


P n c 0 U\ 


(C-18) 


where G = jy/Tn , T and n are given by Eqs. (C-13) and (C-14). The same result can 
be obtained for the relationship between (E 2 ,w,) and( /),»,) . Thus, the transfer matrix 
for the duct element shown is given by 

cosh(rkZ) sinh(r kL)t G 


T = 

1 DE 


Gsinh [VkL) cosh (fkZ.) 


(C-19) 
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(C-20) 


when the wave field in the duct is lossless and isentropic, 

r = / g = l. 


Eq. (C- 1 8) can thus be simplified to 

P 2 

PqCqUi 

or in the p' — Q formulation 

> 2 _ = 

_Qi_ ~ 


cos(kZ) 

ysin(kC) P ] 

ysin(kZ) 

cos (kL) PqC q U x 

cos(kZ) 

j — sin(kZ) r 
S 1 ’ P, 

jPPsm(kL) 

cos(kZ) L2 i. 


0 


(C-21) 


(C-22) 


The determinant of the transfer matrix of Eq. (C-22) is given by 


det 


cos(kZ) j -^-sin(kL) 
$ 

j — sin (kZ) cos(AZ) 


= 1 . 


(C-23) 


Thus the lossless isentropic duct element is a reciprocal system. Furthermore, T n and 
T 22 terms are real, and Jj 2 and T 21 are purely imaginary; thus this system is also 
conservative. 


Clamped piezoelectric backplate with shunt loads 

To develop the transfer matrix representation for the piezoelectric-backplate (or PZT- 
backplate), it is assumed that the maximum dimension of the PZT-backplate is much less 
than the bending wavelength of interest, and that only linear behavior of the PZT- 
backplate need be considered. With these assumptions, an equivalent two-port networks 
including a lossless transformer can be developed for the PZT-backplate, as shown in 
Figure C-2, where P l is the pressure exerting on the PZT-backplate, Q x is the incident 

volume velocity, V and / are the voltage the current at the electrical port, respectively. 
Associated with the two-port networks are two impedances that are measurable properties 
of the system 

= — ^ - — : the blocked electrical impedance (C-24) 

Q = 0 J^^eB 



65 


= R a o + j ( °M aD + : the short-circuit acoustic impedance (C-25) 

joCaD 

where the definition of C eB , R aD , M aD and C aD are listed in Table 1. Furthermore, the 

lossless transformer converts energy between the acoustical and electrical domains and 
obeys the transformer relations 

r=(M2x, (C-26) 



and 

P' = V(j> , (C-27) 

where (f) is the impedance transformation factor. Hence, for two-port network shown in 
Figure C-2, one has 




j(oC eB jcoC eB 
P x =V(f> + Z aD Q x . 

Substituting Eq. (C-28) into (C-29) results in 


P=-^— + Z, 


1 + - 


<t> 


2 \ 


Qv 


(C-28) 

(C-29) 

(C-30) 


eB J 


1 jcoC eB ““I' ja>C t 
where Z aD +<f> 2 / jcoC eB defines the open-circuit (/ = 0) acoustic impedance of the PZT 
backplate. One the other hand, substitution of P ] and V for Q x in Eq. (C-28) from Eq. 
(C-29) results in 


V = - 


/ 


jo>C ei 


j coC eB Z a D 


jcoC eB Z aD + (f> 


+ - 


K 


j(oC eB Z aD + (f> 


(C-31) 


where 


1 

f j<x>C eB Z aD ' 

1 

\ f/UaC.,) 


v. JdCesZaD + ^ , 

c eB 

z, D +f-/(j<»c, B ) 


C 


1 (l-r 2 ) = -^, (C-32) 


eB 


eF 


where C eF is the free electric compliance (/) = 0 ), and k 1 is the coupling factor. Finally, 
from Eqs. (C-28) and (C-29), one has 

A 

<t> <t> 


p _ j^c^Zgp + <t> y ^ aD j 


(C-33) 
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(C-34) 


q = M^!L V 

<!> 



or in the matrix form of 




jcoC eB Z aD + (f> 

1 

i 

N 

1 


~p: 


<t> 

<t> 

~V~ 

_q:_ 


jo)C eB 

1 

I 



</> 

1> _ 



(C-35) 


Consequently, when the electrical boundary condition defined by V/I = -Z eL is given, 
the acoustic impedance, defined by P x /Q x , can thus be obtained for a given PZT- 
backplate as 


p t ( KA + f-)(i / /i)-z„ D ( KA + £_ ) z tl + z, D 
2, jo>e„(vfi)-\ j(oC eB z eL + 1 

The P-p 0 c 0 u formulation of Eq. (C-35) is given by 


Py 


j<X>C eB Z aD + (f> 

<t> 

<t> 

~v~ 

p 0 c 0 U x _ 


jop Q Cq C eB 

Po c o 

i 



A D (f) 

a d </>_ 



where A D is the effective area of the PZT-backplate due to the displacement of the 
clamped PZT-backplate is not unifonn. 
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A 5’ P 5’ U 5 


a,,p,,u 4 


(a) 


APT j APT! 

^ ^2 ’> ± 2 5 ^ 2 


(b) 


L 


P 2 ,t / 2 




C£ 

CvJ 


(c) 


Figure C-l. (a) Schematic of an area contraction, (b) Schematic of an area expansion. 

(c) Schematic of an acoustic duct system. 
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(a) 



(b) 

Figure C-2. (a) Schematic of a piezoelectric backplate with shunt loads, (b) The 

equivalent circuit representation for the piezoelectric backplate with the two- 
port network indicated. 
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Appendix D. Comparison Between TM and LEM 

Equation Section (Next) 

As discussed in Section 4, when the dimensions of the acoustic device are much 
smaller than the wavelength of interest, The device components can thus be lumped into 
idealized discrete circuit elements 16 . In comparison with the transfer matrix 
representation of the acoustic system, the LEM decouples the temporal and spatial 
variables associated the acoustic field, there is thus no spatial variation for a “lumped” 
element. Furthermore, the lumped element is two-pole system, while transfer matrix 
representation is four-pole network. 

On the other hand, the LEM is occasionally no more than a reduced version of the 
transfer matrix representation when the quasi-state assumption is satisfied. For example, 
as shown in Figure D-la, a duct with a sound-soft termination, it is assumed that no 
viscous loss is taken into account, thus from Eq.(C-21), one has 


J g_ = y Wt an« a pi askl<<l 

a ^ s 


(D-l) 


where pJ/S is the lumped acoustic mass of the short tube. One more example, as 
shown in Figure D-lb, the duct is ended by a sound-hard termination, similarly, from Eq. 
(C-21), one has 


A 

a 


-jp 0 c 0 cot kl 
S 


as kl « 1 , 

V 

2 

Po c o 


(D-2) 


where V = SI is the volume of the tube, V/ p Q c 0 is the lumped acoustic compliance C aC . 

However, in most case, more assumptions should be satisfied for the LEM 
consistence with the transfer matrix representation. Shown in Figure D-lc, a duct 
terminated with a complex boundary Z T , from Eq. (C-21), one has 
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P l Z T coskl + j[p 0 c 0 / S) sin kl p Q c 0 Z T (-7 cot kI)/S + (p 0 c 0 /S)~ 

Q ] ( j sin kl ) Z T j ( /? 0 c 0 /S) + cos kl Z T + p 0 c 0 (-j cot kl)/S 

When kl « 1 , Eq. (D-3) reduces to 


(D-3) 


P x _ p 0 c,Z T (-j cot kl)/S + (p 0 cjS ) 2 Z T /(ja>C aC ) + (p 0 c 0 /S ) 2 
Qi Z T +p 0 c 0 (-j cot kl)/S Z T + \/(jcoC aC ) 

where C aC = V/ p 0 cl . Furthennore, when Z T j{ jcoC aC ) » ( p 0 c 0 / S ) 2 , one has 


(D-4) 


Px D z T Kj®c aC ) 

Qx Z T+ \/(jcoC aC y 


which is the LEM of the system. 


(D-5) 
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(c) P 2 /U 2 =Z T 

Figure D-l . (a) Schematic of a duct terminated with sound-soft, (b) Schematic of a duct 
with sound-hard tennination. (c) Schematic of a duct with complex 
tennination 
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Appendix E. Summary On NASA 2DOF Liner Mimic 

Equation Section (Next) 

Introduction 

The efforts have been made to mimic the second layer of the NASA 2DOF liner using 
the compliant diaphragm. The nominal specific acoustic impedance of the second layer 
of the NASA DDOF is shown in Figure E-l. Firstly, the mimic is perfonned using a 
clamped compliant diaphragm to match the target reactance. Proper materials are chosen 
accordingly. Second, a thin PVDF is attached to the compliant diaphragm. In 
comparison with PZT, PVDF is light-weighted and more flexible. A thin PVDF attached 
to the compliant diaphragm will not change it mimic behavior very much. 


Lumped element modeling of a clamped compliant diaphragm 


Unlike a rigid plate, the compliant diaphragm has finite acoustic compliance and 
mass. Next, the lumped element model of a clamped compliant diaphragm is briefly 
introduced. 

Under a uniform pressure P , the deflection of a clamped compliant circular 
diaphragm with radius a and thickness h is given by 


w(r) = w(0) 


/ 



1 - 



V 


J 


(E-l) 


where w(0) is the deflection of the center point of the compliant diaphragm and given by 


w 


( 0 ) = 


Pa 
64 D 


where D , the flexural rigidity, is defined by 


D = 




12(l-v 2 )’ 

where E is the elastic modulus, v is the Poisson’s ratio. 

It is assumed that the diaphragm is under the excitation of a harmonic signal 


(E-2) 


(E-3) 


w(0,t) = w(0)< 


jcot 


(E-4) 
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the physically distributed compliant diaphragm can be lumped into equivalent mass and 
compliance at the center point as long as kd < 1 and the excitation frequency is up to first 
resonant frequency of the diaphragm, where k = <u/c 0 is the wave number and d is the 
diameter of the diaphragm. Firstly, the acoustic mass of the compliant diaphragm is 
determined via computing the kinetic co-energy. The incremental kinetic co-energy 
W . contained in a single ring of the transverse deflection is 


dW k£ , phlnrdr ) | 


jo)w( 0 ) 


Mass 


f 

( r \ 

2 \ 

l- 

— 1 


V 

UJ 

J 


Velocity 


(E-5) 


Thus, the total IT . is 

KE 


a a 

1 

( ( „\ 2 7 

2 

W KE‘=\ dW KE‘=\ n P h 
0 0 

1 

e 

o 

1- ' 

l J 



rdr , 


(E-6) 


thus 


W kf: =^[jcow(0)] 2 ph 


^ na : ^ 


V 5 J 


which must equal the lumped acoustic W 


KE 


L Ma Q 2 =W KE . =^[jo)w(0)J ph 


^ k a 1 \ 


V 5 ) 


where M a is acoustic mass, and Q is the volume velocity 
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( r \ 

27 
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( 2 \ 
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jo)w{ 0) 

1-1 

i 77 1 



2 nr dr = ja>w( 0) 

Kd 
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J 



V 3 J 


The acoustic mass is 


M . = 


9 ph 


2 ‘ 


5 na 

The Eq. (E-9) also indicates that the effective area of the clamped diaphragm is 

, 1 2 

A eff . 


(E-7) 


(E-8) 


(E-9) 


(E-10) 


(E- 11) 
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Then, the specific acoustic mass, M s a , is given by 


M s =M A~ = 

a a ejj 


3ph 


(E-12) 


Next, the acoustic compliance is derived via lumping the total potential energy to the 
center point. Under the pressure P , the incremental potential energy, dW PE is given by 


dW PE = PdAdw(r) = P [2n rdr) d\v{r) . 


Force Deflection 


Moreover, from Eq. (E-l) 


<M')=£$«Mo)= 


Sw(0) 
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UJ 
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dw[ 0) 


and from Eq. (E-2) 


64Dvv( 0) 


(E-l 3) 


(E-l 4) 


(E-l 5) 


So, 


dW PE = 


, , . 128;rZ)w(0) 

= Plnrdrdw ( r ) = — — — 
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rdrdw(O) 


(E-l 6) 


and 
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qw(o] \2%7rDw(0) r 


0 0 


1- 


f-1 

UJ 


2 A 


rdrdw(0) = ^^-\w(0)^ , 


(E-l 7) 


which should equal the lumped acoustic potential energy 

1 1 w2 TJf7 1 64 XD r , x -|2 

V =W PE = 5— w(0) 

2 C„ 2 3 a 1 L V 


(E-l 8) 


where C a is the acoustic compliance, and V is the volume displacement 


V = J w(r)27irdr = J w(0) 
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2\ 

Jw(°) 
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1 „ 1 
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2nrdr = w(0) 


^ na 1 ^ 
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(E-l 9) 


Thus 
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c = 


na 

192 D 


(E-20) 


Then, the specific acoustic compliance is 


c: = 


c 


A 


eff 


64 D 


3(l -v 2 )a 4 

1 


(E-21) 


Eqs. (E-12) and (E-21) give the specific acoustic mass and compliance of the 
compliant diaphragm. Clearly, the specific acoustic mass is the function of the density 
and thickness of the diaphragm, while the specific acoustic compliance relates to the 
radius, thickness, elastic modulus and Poisson’s ratio of the diaphragm. 


Mimic the second layer of NASA 2DOF liner 

To perform the mimic, a curve-fitting method is adopted to obtain set of the specific 

acoustic mass and compliance which satisfy the NASA 2DOF requirements. The curve- 
fitting function is set as 


F = 


j®M s a + 


1 


TAAo) 


(E-22) 


jcoC a j, 

Then, adding ±10% tolerance to the target value of the curve fitting function (see Figure 
E-l), the curving fitting results (thus a mimic region) are obtained as shown in Figure E-2, 
and the range of the acoustic mass and compliance as shown in Figure E-3. 

As mentioned above, the acoustic mass and compliance relate to the geometry and 
material properties of the diaphragm. Then, following the range of the specific acoustic 
mass and compliance, if the geometry parameters of the diaphragm are given, then the 
material properties can be determined. From Eq. (E-12), 

5 M’ 


P 


3 h 


(E-23) 


From Eq. (E-21), 


E = 


3(l-v 7 )a 


1 6/f C 


Thus, if the range of the radius and thickness of the diaphragm are given by 

0.04<?-3 < /? < 2e-3 [m], 
and 


(E-24) 


(E-25) 


76 


2e-3<a<10e-3 [m], (E-26) 

The range of the density of the diaphragm is then obtained as shown in Figure E-4. 
Similarly, the range of the elastic modulus can also be obtained if the range of the radius 
and the thickness of the diaphragm are given. Then, when the range of the geometry 
parameters and material properties are obtained, the next step is to search materials for 
the diaphragm to satisfy the requirement. 

Choice of the materials — Polymer 

Following the guidelines set in section above, one possible polymer material is 
Polytheretherketone (PEEK). The PEEK has high temperature resistance, good 
resistance to wear, dynamic fatigue. The material properties of the PEEK are listed in 
Table E-l. 


Table E-l. Material properties of the PEEK 


Density 

Kg/m 3 

Young’s Modulus 
10 9 N/m 2 

Poisson’s ratio 

1260 - 1320 

3.7-4 

0.34 


Firstly, the diaphragm is not attached with PVDF film. Then, following the range of 
the specific acoustic mass given in last section and Eq.(E-23), the thickness of the PEEK 
film is 

O.le-3 < /z < 0.127e-3 [m], (E-27) 

The thickness of the PEEK diaphragm is chosen as 0.125e-3 m due to the availability of 
the commercial product. Thus, from Eq. (E-23), one has 

M s a = 0.0945 (E-28) 

with the range of the specific acoustic compliance is (shown in Figure E-3) 

2.3e-7<C; <2.9e-7. (E-29) 

The best mimic result happens at 

C s a = 2.46e-7 , (E-30) 

as shown in Figure E-5. From Eq. (E-21), the radius of the PEEK diaphragm is 

a = 10.2e-3 [m] (E-31) 
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Now assuming that the PEEK diaphragm is loaded by a uniformly distributed load of 200 
N/m, the maximum stress is then given by 

fi-L. =Yll il+v)a2 =669eS [Pa] (E_32> 

K»L„=|r^' ,! = 9 ' 99e5 [Pa] (E-33) 

which is much less than the tensile strength of the PEEK (70-100 MPa), 

Next, a thin metalized PVDF is attached to the PEEK diaphragm, as shown in Figure 
E-6. The metalized PVDF film consists of PVDF metalized on both sides with aluminum. 
The thickness of the commercial metalized PVDF is 0.009e-3 m or 0.1 le-3 m. Here, we 
choose the metalized PVDF film with the thickness of 0.009e-3 m. The material 
properties of the metalized PVDF are listed in Table 6. Because the density, Poisson’s 
ratio and Young’s modulus of the PVDF are similar to the PEEK, the thin PVDF film 
attached to the PEEK diaphragm doesn’t change the mimic behavior very much. The 
mimic behavior of the PEEK diaphragm attached with thin PVDF film is shown in Figure 
E-7. The results show that one can mimic the acoustic behavior of the second layer of the 
NASA DDOF liner pretty well. However, further investigation find that the frequency 
shift between open-circuit and short-circuit of the PEEK-PVDF diaphragm is very small 
because the PVDF has relative weak electromechanical transduction in comparison with 
PZT material. Thus, work in the future will focus on the improvement of the frequency 
shift by increasing the thickness of the PVDF or adopting annular PVDF diaphragm 
(ring-like). 

Choice of the materials — Metal 

In comparison with Polymer materials, some light-weighted and flexible metals, such 
as Aluminum/Magnesium alloy, brass alloy and Titanium alloy, are more attractive in 
consideration of the operation temperature and fatigue. Section 6 presents the application 
of Titanium alloy to mimic the 2DOF liner. More metals will be taken into account in 
future. 
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2.4 





Freq.fHz] 

Figure E-2. Curve fitting (mimic region) for target acoustic reactance 
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Figure E-3. Mimic range of the specific acoustic mass and compliance 
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Figure E-4 Range of the density of the diaphragm corresponding to specific acoustic mass range and 

a given thickness range of the diaphragm. 
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Freq. [Hz] 

Figure E-5. Mimic behavior of the PEEK diaphragm 
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Figure E-6. Schematic of composite diaphragm 
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Appendix F. Impedance prediction code - ZKTL 


This MATLAB code is rewritten based on NASA ZKTL Fortran code. The code was 
developed to calculate acoustic impedance for kinds of acoustic liners and the EMHR. 

The code consists of five subroutines listed below: 

• Single channel impedance calculation (linear) — SCIC 

• Single channel impedance calculation (nonlinear) — SCICNL 

• Multi-channel, multi-segment, multi-layer impedance calculation (linear) — 
MCMSML 

• Multi-channel, multi-segment, multi-layer impedance calculation (nonlinear) — 
MCMSMLNL 

• EMHR 

Main menu 

Executing zktl on MATLAB command window results in: 

Main Menu: 

(1) Single Channel - Linear (SCIC) 

(2) Single Channel - Nonlinear (SCICNL) 

(3) Multi-Channel/Segment/Layer - Linear (MCMSML) 

(4) Multi-Channel/Segment/Layer - Nonlinear (MCMSMLNL) 

(5) EMHR - Electromechanical Helmholtz resonator 

Entering selection (1-5) will choose related computing module for acoustic impedance 
calculation for acoustic liner / resonator. Each module has detailed help information. 

SCIC routine 

Choosing (1) on main menu will calculate single channel impedance linearly. In this 
routine, the acoustic impedance of a single channel is independent on the incident sound 
pressure level (SPL) and grazing flow speed. A single channel liner can consist of 
several layers. Each layer consists of one open channel, a perforated plate or foam. 

There are two kinds of geometry for one open channel— tube and split, shown in Figure F- 
1 . The changes of diameter of the open channel from one layer to another layer are 
allowed if the channel is in the form of the tube. For perforated plate, it’s assuming that 
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the plate is infinitely thin. There are two types of input information to define the 
characteristics of the plate. These are: 

• Flow resistance, in this case, the perforated plate only adds an acoustic resistance 
component to the acoustic impedance calculation. 

• Pressure drop across the plate in the complex form, dp{f ) = A(f) + iB(f) where 
A and B are frequency dependent coefficients. 

For a foam layer, the thickness and flow resistance should be included in the input 
parameters. The conclusion for input parameters for SCIC is shown in Figure F-2. 

SCICNL routine 

The option (2) in main menu is to nonlinearly compute the normalized acoustic 
impedance of a single channel liner. SCICNL takes account for the effect of SPL, 
grazing flow and bias flow on the acoustic impedance. In SCICNL routine, a single 
channel also consists of layer(s). Each layer has an open channel or perforated plate or 
both. The foam (layer) is not used in this case. The input information for SCICNL is 
almost the same as the SCIC, except for 

• Grazing flow speed 

• Bias flow speed 

• Input file contains data to define frequency-dependent SPL 

• Coefficients for the definition of changes of impedance across the perforated 
sheet during the nonlinear computing pass. 

MCMSML routine 

Choosing option (3) on main menu will execute MCMSML routine. MCMSML 
routine is used to calculate linearly the acoustic impedance of the liner which consists of 
numbers of channels, segments or layers, shown in Figure F-3. The liner can have a 
group of segments. It is assumed that all segments have the same span- wise width. Each 
segment consists of series of vertical elements. A single element includes open channel 
(s), perforated plate (s) or foam, as same as what described in SCIC routine. The element 
length distribution in each segment can be linear or on Quadratic residue series. The 
input logic for MCMSML is depicted in Figure F-4. 
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MCMSMLNL routine 

The option (4) on main menu is to run MCMSMLNL module to computing acoustic 
impedance of multi-channels/segments/layers liner using nonlinear acoustic method. The 
liner geometry is the same as what is depicted in Figure F-4, except for that the foam 
(layer) is not used in this case. In addition to what are included in the input information 
for MSMSML, some extra input parameters are needed to define MSMSMLNL routine: 

• Bias flow velocity 

• Grazing flow velocity 

• Input file which contains data to define frequency-dependent SPL 

• Coefficients for the definition of changes of impedance across the perforated 
sheet during the nonlinear computing pass. 

EMHR routine 

Choosing option (5) on main menu will execute EMFIR routine. The EMHR routine 
uses LEM or TM to predict the acoustic impedance of the EMHR. The input parameters 
include: 

• Dimensions of the EMHR 

• Material properties of the piezoelectric composite backplate 
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Figure F-l: Channel geometry for SCIC/SCICNL 



Figure F-2: Input parameters for SCIC 
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(?) : Open channel length of 
segment 1 at x_min 

(2) : Open channel length of 
segment 1 at x_max 

d) : Open channel length of 
segment 2 at x_min 

0) : Open channel length of 
segment 2 at x_max 
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Figure F-3: Liner geometry for MCMSML/MCMSMLNL 
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Figure F-4: Input logic for MCMSML 
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